Silicon surface with giant spin-splitting 
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Abstract 

We demonstrate the induction of a giant Rashba-type spin-splitting on a semiconducting substrate by 
means of a Bi trimer adlayer on a Si(lll) wafer. The in-plane inversion symmetry is broken so that the 
in-plane potential gradient induces a giant spin-splitting with a Rashba energy of about 140 meV, which is 
more than an order of magnitude larger than what has previously been reported for any semiconductor het- 
erostructure. The separation of the electronic states is larger than their lifetime broadening, which has been 
directly observed with angular resolved photoemission spectroscopy. The experimental results are con- 
firmed by relativistic first-principles calculations. We envision important implications for basic phenomena 
as well as for the semiconductor based technology. 
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Exploiting the electron spin for information processing is one of the leading goals in the rapidly 
growing field of spintronics. At its heart lies the Rashba-Bychkov (RB) type spin- splitting, where 
the spin-orbit interaction lifts the spin degeneracy in a symmetry broken environment Q2Q . Man 
device proposals make use of this concept |2> siSSl with some interesting proofs of principle 
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. The materials of choice are semiconductor heterostructures, albeit the size of the spin- splitting 
is typically very small. A large spin- splitting is desirable as it would, for example, decrease the 
precession time of the spin in a spin transistor yD so that it is smaller than the spin relaxation time. 
Furthermore, a separation of the spin-split states beyond their lifetime broadening is an important 



The 



criterion for distinguishing between the intrinsic and extrinsic spin Hall effect 
different influences on the intrinsic spin Hall conductivity, such as disorder and elastic/inelastic 
lifetime, are still under debate [Jil llL 



_ 3.ecently f a giant spin- splitting has been demonstrated for noble metal based surface alloys 
OH ll4j,li5D, where heavy elements with a strong atomic spin-orbit coupling are incorporated 
into the surface. These systems, however, are not suitable for the field of spintronics because of 
the presence of spin-degenerate bands at the Fermi level originating from the metallic substrate. 
One possible alternative is to grow thin films with spin-split bands onto a semiconducting substrate 
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18j]. However, due to confinement effects a multitude of quantum well states arise, which 



potentially influence the transport properties of the system. It is, therefore, desirable to transfer 
the concept of the giant spin- splitting directly onto a semiconductor surface. 

Here we show that a monolayer of Bi trimers on a Si(lll) surface forms a two-dimensional 
(2D) electronic structure with a giant spin- splitting much larger than what has been observed 
so far at the interfaces of semiconductor heterostructures. The effect can be traced to a strong 
contribution of an in-plane potential gradient due to an inherent structural inversion asymmetry 



(RB model). Whi 
experimentally [1, 



e the structure of this system has been studied both theoretically as well as 
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2 ill , the electronic structure, in 



particular a possible spin-splitting of the 



24fl . We demonstrate unequivocally that 



electronic states, has remained a controversial issue [23, 
Bi induces a giant spin-splitting at the silicon surface. Furthermore, the spin- splitting is observed 
to be larger than the lifetime broadening, so that the Bi/Si(l 11) system is a prime candidate for 
spintronics applications or studying the intrinsic spin Hall effect. In addition, the silicon substrate 
allows for excellent compatibility with existing silicon-based semiconductor electronics. 

A single layer of Bi on Si(lll) grows in a monomer as well as a trimer configuration, both 
of which show a (\/3 x v / 3)R30° reconstruction (iLllSLIZlD- A structural model is shown in Fig. 
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FIG. 1: (color online) Structural model of the two (V3 x V3)R30° phases of Bi/Si(lll): (a) monomer 
phase (b) trimer phase. The thin black lines indicate mirror planes of the Bi adlayer. The thicker black lines 
indicate the (\/3 x \/3)R30 o unit cell. The smaller the spheres, the further away they are from the surface. 

CD for the monomer phase (a) and the trimer phase (b). Both the monomers and the trimers are 
centered on top of second layer Si atoms (T 4 lattice sites) forming a symmorphic space group 
based on the point group 3m. The Si substrate breaks the in-plane inversion symmetry for both 
the monomer and the trimer phase. Looking at the isolated Bi adlayer alone, the trimer formation 
introduces a reduction of the symmetry because the mirror plane a v2 is missing. The mirror plane 
cr„i holds for both the monomer and the trimer phase as well as for the combination of adlayer and 
Si substrate. From these simple symmetry considerations we conclude that the Bi trimer phase is 
the least symmetric structure and, hence, should lead to the bigger spin-splitting. We, therefore, 
only consider the trimer phase. Its preparation was verified with quantitative low-energy electron 
diffraction measurements Il22ll . 

The experimental band structure measured with angular resolved photoemission spectroscopy 
(ARPES) along the two high symmetry directions of the surface Brillouin zone (SBZ) T M and 
T KM is displayed in Fig.|2](a) and (b), respectively D22Q. The intense feature near T at an energy 
of about —2.3 eV can be attributed to the silicon bulk. The other features (SI, S2, S3) in Fig. (2a) 
originate from the 2D electronic structure of the surface. SI is most intense at the M-point at an 
initial state energy of about —1.3 eV. This band splits in two components when moving away from 
the high symmetry point M, which is a strong indication of a RB-type spin- splitting. S2 is located 
at about — 2.3eV at the second T-point and disperses upwards towards the M-points. The third 
state S3 shows the highest intensity at the second M-point at an energy of about — 2.5eV. This 
band moves downwards in energy towards the second T-point. The bandwidth of S3 is smaller 
than the one for S2. These three 2D states are also visible along the T K M-direction as shown in 
Fig. |2](b). SI appears as a parabolic band with negative effective mass with a band maximum 
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FIG. 2: The two panels show angle-resolved ultra violet photoemission spectroscopy data of the trimer 
phase of bismuth on Si(l 11) along the two high symmetry directions T M The photoemission intensity is on 
a linear scale with black and white corresponding to highest and lowest intensity, respectively. The energy 
scale is set to zero at the Fermi level, (a) and TKM (b). Light an dark areas correspond to low and high 
intensities, respectively. A splitting of the two-dimensional state into two bands around the M point along 
the T M-direction at an initial state energy of about — 1.3eV is clearly visible in panel (a). We attribute 
this splitting to the Rashba-Bychkov effect with a momentum offset ko = 0.126 A -1 and a Rashba energy 
E R = 140 me V. 



located at about — 1.3 eV at the M-point. Along the TKM direction no splitting of this band has 
been observed. S2 is located around — 1.8eV at M, but only with a very weak intensity. The 
most intense feature along the T K M-direction is the S3 2D state with a band minimum at about 
—2.5 eV at the M-point and an upwards dispersion towards the neighboring K-points. 

A possible spin- splitting in the Bi/Si(l 11) system is an unresolved issue in the literature. While 
Kinoshita et al. [|23|l consider a splitting in the three 2D states related to a strong spin-orbit interac- 
tion of the Bi atoms, it has been dismissed by Kim et al. 112411 . In the following, we will show from 
the experimental data as well as spin-resolved band structure calculations that the band structure 
shows a giant spin- splitting of the electronic states due to the RB effect. 
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FIG. 3: Experimental band structure of Bi/Si(lll) near the M-point. The measurements along TMT (a) 
and KM K (b) show the anisotropic topology of the spin-split bands. 



Spin-degeneracy is a consequence of both time reversal and spatial inversion symmetry. 
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latter is broken spin-degeneracy can be lifted by the spin-orbit interaction (RB model) [22]. In 
addition, if the inversion symmetry is also broken in the plane of the two-dimensional electron 



gas, the contribution from an in-plane gradient can strongly enhance the spin- splitting II 1 311 . The 
characteristic parameters quantifying the strength of the splitting are the momentum offset k , the 
coupling constant in the Hamiltonian an (Rashba parameter), as well as the Rashba energy Er. 
They are related by Er = h 2 k 2 l /2m* and k = nfau/h 2 . Here m* is the effective mass. 

A close up of the band structure near the M-point is shown in Fig. [3] The bands along T M (Fig. 
Oa)) near — 1.2 eV clearly show the characteristic dispersion of a RB type spin- splitting with the 
band crossing at the M-point and the shift of the maxima away from it. From the data we extract 
the momentum offset k = 0.126 A -1 , an effective mass of m* = 0.7 m e (m e free electron mass) 
as well as the Rashba energy Er = 140 me V. From these values we can calculate the Rashba 
parameter a#=1.37eVA. The spin- splitting is well resolved in the data. The average line width 
for the spin-split states at the band maximum (k x = — 0.126 A -1 ) is 195 me V, which accounts 
for intrinsic lifetime as well as interactions and scattering. The separation of the states is about 
220 meV. 

The spin- splitting at the M-point in Fig.[3£a) is strongest along the T M-direction. Along the 
K M K-direction in Fig. [3£b) the spin- splitting at the M-point is much weaker and cannot be re- 
solved in the experiment. This peculiar band topology can be related to the symmetry properties 
of the M-point. As the M-point is located on the border of the first SBZ it has a reduced symmetry 
as compared to the T-point. Despite the symmetry breaking of the Bi trimers, the mirror symme- 
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FIG. 4: (color online) Theoretical band structure calculations for the trimer phase of bismuth on sili- 
con^ 11). Panel (a) and (b) show the calculated dispersion along TM and TKM, respectively. Blue and 
red correspond two opposite spin-polarizations. The calculated spectra reproduce the main features of the 
measured band structure, especially the spin-splitting of the bands around the M-point along T M. 

try a v i (see Fig. [B holds so that for the dispersion along the K M K-direction the spin- splitting is 
greatly reduced, i. e. it can not be observed in the data. 

To support our interpretation of the observed spin-splitting, we conducted spin-resolved first 
principles band structure calculations, which were performed in close analogy to our previous 
calculations on the RB effect H 1 3M . The surface geometry of the trimer structure is determined from 
first principles using the VIENNA AB-INITIO SIMULATION PACKAGE (VASP) which provides 
precise total energies and forces |25|]. The Bi trimers (milkstool structure) are relaxed outward by 
13 % from the ideal positions (100% corresponds to the Si bulk interlayer distances, lattice constant 
5.403 A). The subsurface relaxations are small (< 0.5%) and neglected in the Korringa-Kohn- 
Rostoker (KKR) calculations. The in-plane displacement of the Bi trimer atoms S is 0.3, with 
(5 = indicating Bi on-top of first layer Si atoms and 5 = 1 coinciding Bi-trimer atoms on T 4 sites. 
The subsequent KKR and relativistic layer-KKR calculations use the structural data from VASP 
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as input. The spectral density n±(E, feu) is obtained from the imaginary part of the site-dependent 
Green function. Resolved with respect to spin orientation (index ±) and angular momentum, it 
allows a detailed analysis of the electronic structure. The difference n + (E, fey ) — n_ (E, fey ) reveals 
the characteristic spin splitting of RB-split bands. 

The results of the band structure calculations are shown in Fig. |4] for the T M-direction in (a) 
and for the T K-direction in (b). The intensity scale shows the total spectral density (n + (E, fen) + 
ri-(E, feii)) of the states multiplied by the sign of the spin-polarization sgn(n + (E, fen)— ri-(E, fey)) 
perpendicular to the high symmetry line, i. e. blue and red colors correspond to opposite spin- 
polarizations. The calculations reproduce all the main features of the measured band structure. 
In particular, the splitting of the SI band around the M-point along the T M-direction is well 
documented. As can be seen in figure[4]the two branches of the split S 1 band clearly show opposite 
spin-polarization, i. e. a giant spin- splitting in the electronic structure of Bi/Si(l 11). 

The spin- splitting is strongly anisotropic around M. The peculiar band topology, which was 
observed in the experiment is clearly reproduced in the calculations. This can again be attributed 
to the lower symmetry of wave vectors k\\ within (r M) or perpendicular (r K M) to a mirror plane 
of the system. It is conceived that this feature results from the 'trimerization' of the three Bi sites 
in the 2D unit cell; calculations with a reduced 5 (i. e. larger distance between Bi trimer atoms) 
indicate an even smaller splitting along T K M. Furthermore, the calculations show that about 83% 
of the spin-split states at the M-point are localized in the Bi adlayer and about 16% in the first Si 
layer. One can thus speculate that the spin- splitting is particularly influenced by the Bi adlayer and 
that trimerization symmetry breaking increases the effect of the in-plane potential gradient. 

The giant spin- splitting in the Bi/Si(l 1 1) trimer system has a similar origin as in the Bi/Ag(l 1 1) 
surface alloy: An inversion symmetry breaking in the plane of the surface leads to a strong contri- 
bution from an in-plane potential gradient, which substantially enhances the spin- splitting. In both 
systems the threefold symmetry of the underlying substrate breaks the in-plane inversion symme- 
try. However, considering only the topmost layer, the trimer formation in Bi/Si(l 11) also leads to a 
breaking of the in-plane inversion symmetry (see Fig. [I]), which is not the case for the Bi/Ag(l 1 1) 
surface alloy. 

Comparing the spin- splitting of the Bi/Si(lll) electronic structure to semiconductor het- 
erostructures, we find that in the latter the spin- splitting is substantially smaller. For example, 
for an inverted InGaAs/InAlAs heterostructure a Rashba constant of cxr = 0.07 eVA has been 



measured Q26Q . With an effective mass of m* = 0.05 m e , a Rashba energy of Er = 16 fieM can 
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be calculated. For HgTe quantum wells a Rashba constant olr = 0.45 eVA has been found [27]. 
However, here the spin- splitting has been identified to be proportional to fcn instead of a linear 
dependence [28]. For the Bi/Si(l 11) system, the Rashba energy E R = 140 meV as well as the 
Rashba parameter = 1.37 eVA are much bigger. From the momentum offset k = 0.126 A -1 
we can calculate that a phase shift of the spin precession angle A9 = n can be obtained after a 
length L = A9/2k of only 1.3 nm. In the InGaAs/InAlAs hetero structure a length of 400nm 
has been estimated. While these figures show the excellent potential of the Bi/Si(l 11) system, 
additional measurements giving insight into the transport properties, such as Shubnikov-de Haas 
oscillations, are necessary to further elaborate the suitability of this system for spintronics appli- 
cations. Corresponding experiments are in progress. 

We have shown that the trimer phase of Bi on Si(l 11) shows a giant spin- splitting. The experi- 
mental results reveal the characteristic band dispersion of a RB-type spin-splitting with a peculiar 
band topology at the M point. They are confirmed by first principles band structure calculations. 
The splitting is caused by the spin-orbit interaction induced RB effect in combination with a strong 
contribution from the in-plane gradient due to the reduced symmetry of the trimer structure and the 
substrate. Furthermore, this spin- splitting is of the same order of magnitude as the one reported 
for Bi/Ag(lll) and orders of magnitude lager than a typical spin-splitting reported for semicon- 
ductor heterostructures. In this way, we have transferred the concept of giant spin- splitting onto 
a semiconducting substrate. This gives excellent perspectives for the use of this concept in the 
field of spintronics. In particular, the silicon substrate makes this system compatible with existing 
semiconductor technology. On the fundamental side such systems are interesting for, e. g., the 
spin Hall effect. Since the energy separation of the spin-split states (220 meV) is larger than the 
lifetime broadening (195 meV), it may be easier to distinguish the extrinsic and intrinsic spin Hall 
effects. 
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Supporting Information 
Sample preparation 

The sample preparation as well as the measurements were conducted in ultra high vacuum 
(UHV) with a base pressure of 2 x lCT 10 mbar. The n-doped Si(lll) substrate was annealed at 
1100°C by direct current heating for 10 minutes and cooled down slowly to 800°C over a time 
interval of 10 minutes until a sharp (7x7) low energy electron diffraction (LEED) pattern was 
observed. One monolayer (ML) of Bismuth was deposited at a substrate temperature of 470°C 
using a commercial electron beam evaporator to form an adlayer of trimers. The substrate temper- 
ature was measured with an optical pyrometer. After Bi-deposition LEED measurements showed 
a(V3x v/3)R30 o reconstruction. 

Quantitative LEED measurements 

As both the monomer phase at 1/3 ML Bi coverage and the trimer phase at 1 ML Bi coverage 
show the same (a/3 x v / 3)R30° reconstruction we used quantitative LEED measurements to dis- 
tinguish between the two phases. We measured the integrated intensity of the (10) and (01) spots 
as a function of electron energy and compared them to calculations done by Wan et al. [1]. The 
measured data was averaged over three equivalent spots and smoothed. The result is shown in 
figure [51 The agreement between measured and calculated spectra is quite convincing, allowing 
an unambiguous identification and preparation of the trimer phase. 

STM measurements 

In figure [6] two (10 x 10)nm 2 topographic STM images are shown. In (a) the monomer phase is 
of the Bi/Si(l 1 1) structure is shown at a bias voltage of —1.55 V and a tunneling current of 0.2 nA 
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FIG. 5: The bismuth coverage was determined via quantitative LEED measurements. The figures show the 
integrated intensity as a function of electron energy for the (10) and (01) spots of the monomer and trimer 
phase, respectively. 

(b)§i/Si"llT)trimers 

FIG. 6: Topographic lOxlOnm 2 STM images of the (-y/3 x ^3)R30 o phase of Bi/Si(lll). The monomer 
phase and the trimer phase are shown in (a) and (b), respectively. 

is shown. In (b) the trimer phase is shown at a bias voltage of —0.95 V. The tunneling current was 
0.2 nA as well. 

Band structure measurements 

Angle-resolved photoemission spectroscopy (ARPES) measurements were done with a hemi- 
spherical SPECS HSA3500 electron analyzer with an energy resolution of ~10meV. The images 
were recorded with a step size of 1° in angular direction. We used monochromatized He I radiation 
with an energy of 21.2 eV for the ARPES measurements. The sample was kept at 90 K during the 
measurements. We measured the band structure along the two high symmetry directions TM and 
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FIG. 7: Sketch of the surface Brillouin zone (SBZ): the blue and black hexagons show the (1 x 1) SBZ of 
the Si substrate and the (y3 x \/3)R30 o SBZ of the trimer phase of Bi/Si(l 1 1), respectively. The red arrows 
indicate the two high symmetry directions TM and T KM with respect to the (\/3 x \/3)R30 o SBZ. 

T KM of the surface Brillouin zone (SBZ) associated with the (v^3 x V3)R30° reconstruction (see 
Fig. ED. 

The zero for the initial state energy in the measured data was set to the Fermi level. We chose an 
n-doped Si substrate because the position of the Fermi level is close to the conduction band. This 
allows states within the bulk band gap to be easily accessible by photoemission spectroscopy. From 
the recorded data we subtracted a Shirley background and normalized every energy distribution 
curve (EDC) of the image to the same integrated intensity. 



Rashba-Bychkov model 

Spin-degeneracy is a consequence of both time reversal and spatial inversion symmetry. If the 
latter is broken spin-degeneracy can be lifted by the spin-orbit interaction. For a free electron 
gas in two dimensions the spin-orbit coupling Hamiltonian in the Rashba-Bychkov (RB) model is 
given by the following equation Q2Q : 



H S o = a R cr(k n x e z ) 

where the coupling constant ajt is called Rashba parameter, a are the Pauli matrices, ku is the 
in-plane momentum and e z is a unit vector normal to the plane of the two dimensional electron 
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gas (2DEG). The resulting energy dispersion is: 

with the effective mass m*, the offset in parallel momentum k and the band externum at E . 
The RB Hamiltonian lifts the spin-degeneracy and introduces an offset in parallel momentum of 
the nearly free electron parabola. At high symmetry points the bands are degenerate due to time 
inversion symmetry. As a result the band dispersion of a RB split band has a very characteristic 
shape. The Rashba parameter an in the most straight forward formulation of the RB Hamiltonian 
is proportional to the potential gradient perpendicular to the plane of the 2DEG. An additional 
contribution from an in-plane gradient can strongly enhance the spin- splitting. 

From the photoemission data the Rashba energy E R (the energy difference between the crossing 
point of the two parabolas and the band externum), the offset of the two parabolas in parallel 
momentum k and the effective mass m* can be determined. These parameters are related to an 
effective Rashba-parameter an via 

En = ^ = ~ 2 a R k (1) 

This effective Rashba parameter cxr contains all the different contributions to the spin-splitting, 
i. e. the potential gradient in the plane of as well as perpendicular to the surface and the atomic 
spin-orbit coupling from the atoms involved. The Rashba parameter cur determined from equation 
(OQ) should be understood as a characteristic parameter for comparison with other systems. 
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